Pairing dynamics in strongly correlated superconductivity 
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Confirmation of tfie phononic origin of Cooper pair formation in superconductors came with tfie 
demonstration tfiat the interaction was retarded and that the corresponding energy scales were 
associated with phonons. Using cellular dynamical mean-field theory for the two-dimensional Hub- 
bard model, we identify such retardation effects in d-wave pairing and associate the corresponding 
energy scales with short-range spin fluctuations. We find which frequencies are relevant for pairing 
as a function of interaction strength and doping and show that the disappearance of supercon- 
ductivity on the overdoped side coincides with the disappearance of the low energy feature in the 
antiferromagnetic fluctuations, as observed in neutron scattering experiments. 



In ordinary superconductors, the origin of attraction 
between electrons, the "pairing glue" , manifests itself in 
observable quantities. Indeed, the characteristic frequen- 
cies of phonons appear directly in the frequency depen- 
dence of the gap function, which in turn enters observ- 
ables such as the single-particle density of states or the in- 
frared conductivity. Migdal-Eliashberg theory P, [2] has 
been extremely successful to extract from these experi- 
ments the spectral function of the phonons that provide 
the glue. 

High-temperature superconductors, heavy-fermion 
and layered organic superconductors all have phase di- 
agrams where non s-wave superconducting order param- 
eters lie in close proximity to antiferromagnetic phases. 
In the case of high-temperature superconductors, much 
effort has been devoted to find out whether antiferromag- 
netic fluctuations could be the pairing glue [1, 0, S, S 0] ■ 
Even though its assumptions are not generally valid in 
that case, Eliashberg theory has been used to extract 
the amplitude and frequency dependence of a spectral 
function that is found to be similar to that for antiferro- 
magnetic fluctuations directly measured by neutron scat- 
tering @,i,[3,|Ti|. 

But understanding the origin of pairing in high- 
temperature superconductors requires an approach that 
does not rely on the assumptions entering Eliashberg 
theory and that takes into account Mott insulating be- 
havior. This seems to rule out theories that are based 
purely on early weak-coupling ideas of boson exchange 
[1^ 13, In fact, Anderson [l^ has argued that the 
appropriate starting point consistent with Mott physics 
is the strong-coupling version of the Hubbard model, or 
the t — J model, and that in this case interactions are 
instantaneous, as suggested by mean- field factorization 
0. 

In this paper, we show, for the Hubbard model, that 
spectral features of the imaginary part of the anomalous 
(off-diagonal) self-energy do correspond to those of the 
spectral function for short-range spin fluctuations and 
that the energy scales relevant for pairing are also those 
of spin fluctuations. 



The Hubbard model Hamiltonian is given by 



H 



E 



(1) 



where tij and U correspond to the hopping matrix and 
the onsite screened Coulomb repulsion respectively with 
c'l'l the destruction (creation) operators for an electron at 

site i with spin a and — c| ^Ci^a the number operator. 
The theoretical method that has been most successful 
to date to treat the Mott transition starting from the 
one-band Hubbard model is dynamical mean-field theory 
(DMFT) 17]. 



Cluster generalizations of DMFT [18|, [IJ, l20|, l2l| are 
necessary to study problems in two dimensions where cor- 
relations beyond single site must be taken into account to 
study, for example, d-wave superconductivity. They lead 
to phase diagrams that have the same features as those 
observed experimentally for both electron- and hole- 
doped high-temperature superconductors 23, IE M, 2^ 
and for organic conductors. In addition, observable quan- 
tities such as the density of states, [24] the ARPES spec- 
trum 2^, 23] and the optical conductivity |23] have the 
experimentally observed behavior. The method that we 
use, C-DMFT with exact diagonalization at T = 0, is de- 
scribed in Refs. and in Ref. [26|. We stress that 
it does not involve any Eliashberg-like approximation. 

In C-DMFT, antiferromagnetism and d-wave super- 
conductivity coexist over part of the phase diagram. This 
is seen in stoichiometric cuprates with intrinsically doped 
planes P] and in a few other cases, [23, [2!! but does not 
appear to be a completely generic property of the phase 
diagram. Our work is thus restricted to showing that spin 
fiuctuations are relevant for high-temperature supercon- 
ductivity all the way to the overdoped regime, leaving 
open the possibility that additional types of fiuctuations 
may either contribute to or hinder superconductivity in 
the underdoped phase. 

The correspondence between the imaginary part of the 
anomalous self-energy, S"„, and the imaginary part of 
the local spin susceptibility, x'\ is seen in Fig. [H We 
take band parameters appropriate for La2_a;Sr3;Cu04, 
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FIG. 1: (Color online) (a) Imaginary part of the anomalous 
self-energy Im Ean = S"„ at the Fermi wave vector nearest 
to the antinodal point, for various dopings. In (b) Imaginary 
part of the local spin susceptibility Imx = x"- Black dots 
in (a) and (b) identify peaks. The position of the peaks of 
E"„ in (a) are reported as magenta dots in (b) at the same 
height as the corresponding x to illustrate the correspon- 
dence between the main peaks of the two functions. The 
frequency splitting between the peaks decreases with doping, 
like the single particle gap. The red curves are for the normal 
state. The lower frequency peak present in the supercon- 
ducting state disappears and the next peak moves to higher 
frequency with doping. In all the figures Lorentzian broaden- 
ing is 0.125t, U = 8t, i' = -0.3f, t" = -O.OSf, f = = 1, 
for La2-2,Sra;Cu04. 

namely t' = —O.lTt for nearest-neighbor and t" = O.OSt 
for next-nearest-neighbor hopping. CDMFT with U = 8t 
then leads to superconductivity in the doping range ob- 
served experimentally [23| . The anomalous self-energy 
E"„ is defined as minus the off-diagonal part of the in- 
verse Green function in Nambu space. Numerical results 
are presented in energy units where t — 1. For all differ- 
ent dopings, the positions of the first two peaks in the 
spin fluctuations (black dots on middle panel) are just 
shifted down with respect to the corresponding peaks in 
S"„ (black dots on top panel). 

In Eliashberg theory for the electron-phonon interac- 
tion, the first two peaks in the phonon density of states 
are shifted down with respect to those in S"„ by the 



BCS gap Similarly, the down shift of peaks in x" 
seen in Fig. [TId increases as we underdope, like the single- 
particle gap. For U — 12t and realistic band structure for 
YBa2Cu307_2; the shift is very weakly doping dependent 

In Migdal-Eliashberg theory, the real part of the self 
energy EJj„ times the quasiparticle renormalization fac- 
tor is the gap function. We find that this function, has 
no static (frequency independent) contribution, 12611 con- 
trary to what was found in the t — J model. [7|, |24| . 

To identify the energy scales relevant for the pairs, we 
introduce the function 

IgH^~ T—lmF^^iu;'). (2) 
Jo 

is the retarded Gork'ov function defined in imaginary 
time by Fij = — (Tci|(r)cjx (0)) with i and j nearest- 
neighbors. The infinite frequency limit of Iq (w) is equal 
to (cj|Cj|) which in turn is proportional to the T — 
d-wave order parameter (it changes sign under 7r/2 ro- 
tation). It was shown in Ref. [24| that (cifCji) scales 
like Tc- For all these reasons, Iq (w) is useful to esti- 
mate the frequencies relevant for binding. Its meaning is 
illustrated by the d-wave BCS result in Fig. [2^). The 
function Iq {lu) rises monotonically until it reaches the 
sharp BCS cutoff frequency lUc above which no binding 
occurs. Ig {(^) extracted from the Eliashberg calculation 
[sot fo'^ Isad is also displayed in Fig. [5^). The maximum 
is reached at a frequency just above the largest phonon 
frequency. 

Ig (w) is plotted in Fig. [^b for underdoping S = 0.4, 
optimal doping S = 0.16 and overdoping S = 0.26. The 
asymptotic large frequency value of Ig (w) indicated by 
horizontal lines gives the order parameter that, as a func- 
tion of doping, has the dome shape dependence [lH . The 
functions Ig (w) cross their respective asymptotic values 
at progressively lower frequencies as doping increases. 
The spin fluctuations that dominate at the lower frequen- 
cies come from wave vectors around (tt, tt) , as illustrated 
in Fig. for an underdoped case. The maximum of 
Ig (w) is more pronounced in the underdoped regime. 
The form of Ig (w) in the overdoped regime is closer to 
the BCS limit with just a sharp cutoff. Our calculations 
are less precise at high frequencies, but nevertheless they 
suggest that, in all cases, Ig (w) undershoots very slightly 
its asymptotic value and then recovers at frequencies that 
are of order U/2 where the upper Hubbard band opens 
new scattering channels [3]. This has no analog in ordi- 
nary superconductors. 

In Fig. [3] we focus on the low- frequency behavior. On 
the top panel, Ig {i-^) crosses its asymptotic value for 
the first time near its maximum. This crossing point 
shown by vertical lines follows the first peak in the corre- 
sponding x" in the bottom panel. By studying the cases 
[/ = 8, 12, 16 we have verified that these features scale 
with J. Clearly, if we wished to design an approximate 
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FIG. 2: (Color online) (a) The dashed green line is la (i^) 
for a d-wave BCS superconductor with a cutoff at Uc — 0.5. 
The maeenta line is obtained from Eliashberg theory for Pb 
in Ref. |30|- Frequencies in that case are measured in units 
of the transverse phonon frequency. The two glitches before 
the maximum correspond to the transverse and longitudinal 
peaks in the phonon density of states. The scale of the verti- 
cal axis is arbitrary, (b) Iq (i^) calculated for various dopings. 
The horizontal lines for the asymptotes mark the value of the 
order parameter, (c) The three independent Fourier compo- 
nents of on a 2 X 2 plaquette for an underdoped case. The 
(tt, tt) component dominates at low frequencies. 



mean-field theory [sT] for this problem that would play 
a role analogous BSC theory as an approximation of the 
Migdal-Elisahberg theory, we would use a frequency cut- 
off of order J. 

Let us now discuss how the properties of the spin 
fluctuations x" i'^) correlate with those of the d-wave 
superconducting state for La2-2;Sra;Cu04. In Fig. [TJ), 
one sees that in the underdoped regime the low fre- 
quency peak is the most prominent feature. Optimal 
doping ((5 ^ 0.16) is reached when the intensity of the 
low frequency peak becomes comparable to the next one 
at higher frequency. More importantly one sees that, 
around doping S = 0.26 , superconductivity disappears 
with the low frequency peak in x" i'^) below lu about 
J/2 = 0.25t. That low frequency peak is the one involved 
in the dynamics of the pairs as shown by the intersection 
of the Ig (w) function with its asymptotic value in Fig. 
^p. The leftover peak in the normal state, indicated in 
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FIG. 3: (Color online) (a) Low frequency part of Iq (i^) 
for three dopings, underdoped, optimally doped and over- 
doped. The vertical lines indicate the first intersection with 
the asymptotic value. The plot of x" in (b) shows that 
the first peak occurs in the frequency interval where Ig (<^) 
reaches its maximum. These are the main frequencies that 
give rise to pair binding. 



red in Fig. [TJd, moves to higher frequencies as we dope. 
These features are those found in neutron scattering ex- 
periments, [s^lsst 

The large tails and monotonic decrease of the weight 
of the low frequency peak in x" dope are con- 

sistent with the "glue function" extracted from recent 



optical conductivity experiments. 11| The position of 
the low frequency peak near 0.2t at optimal doping is 
consistent with the experimental value [llj if we take 
t — 250 meV. One should recall that Fig. [T)d for the lo- 
cal spin spectral weight x" gives information integrated 
in wave vector so the properties of the "neutron reso- 
nance" located at (tt, tt) have to be found by a different 
approach. In a recent calculation with a related cluster 
method, fs?] it has been found that the peak located at 
(tt, tt) in the infinite lattice decreases with frequency in 
the underdoped regime. Given the small weight of this 
"neutron resonance" , this does not contradict the fact 
that x", whether local or averaged over one-quarter of 
the Brillouin zone near (tt, ttL has the opposite doping 
dependence m, M, M^M, m The magnetic resonance 
itself has small weight [38[ . 

We stress that, despite the similarities, the results ob- 
tained in this paper are not identical to those that are 
obtained in ordinary superconductors. In particular, in 
the underdoped regime the spin fluctuations are strongly 
pair-breaking. The pair-breaking effect of the pseudogap 
can be seen from the fact that in the normal state, the 
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pairing susceptibility decreases as one approaches half- 
filling when vertex corrections are neglected [2^ . Follow- 
ing the suggestion of Ref. [s^ , we also checked whether 
the pseudogap is pair breaking by computing the ratio 
Tian (ii^n) / (1 — S„ {iujn) / i^n) Sit the antinodal Fermi sur- 
face crossing as a function of Matsubara frequency. The 
ratio is far from constant, in agreement with the exis- 
tence of strong pair-breaking effects in the pseudogap. 
From a diagrammatic point of view, spin fluctuations 
both scatter electrons (self-energy), decreasing the den- 
sity of states at the Fermi surface, and provide the glue 
(vertex). This can lead to the dome shape of the transi- 
tion temperature as a function of doping 40] . In the un- 
derdoped regime, the pair-breaking effect wins over the 
glue provided by the vertex, whereas in the overdoped 
regime the vertex dominates. 



In conclusion, we have found that the imaginary part 
of the anomalous self-energy has a structure that is cor- 
related with the spectrum of spin fluctuations. This cor- 
relation is similar to the one found with the phonon spec- 
trum in the Migdal-Eliashberg theory of ordinary super- 
conductors. This suggests the importance of spin-one ex- 
citations in pair formation. Our approach also allows for 
mutual feedback between spin fluctuations and pairing 
4l( . The frequencies most relevant for the pair dynamics 
scale with the Heisenberg exchange J. Superconductivity 
disappears for sufficient overdoping when the first peak 
in x" below frequencies about J/2 becomes negligible. 
There are, however, major differences with ordinary su- 
perconductors. Even though the anomalous self-energy 
increases as we approach half-filling, the order parame- 
ter decreases because of large self-energy effects in the 
normal part of the propagator. These come from Mott 
Physics at half-filling. The magnetic fluctuations that we 
find have the same doping and energy dependence as that 



found in optical, [9|. llO 
experiments 



32, m, 35 



111 t unneling [3, |4^ and neutron 
36l. 37 1 . This work leaves open 



the possibility that in the underdoped regime there exists 
other instabilities that compete with antiferromagnetism 
and d-wave superconductivity. 
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METHODS, SUPPLEMENT 



method) Q : 



In Cellular Dynamical Mean-Field Theory (C-DMFT) 
a cluster is embedded in a bath of non-interacting elec- 
trons that simulates the effect of the rest of the infinite 
lattice by injecting and removing electrons on the clus- 
ter with the appropriate single-particle propagator. The 
bath is determined self-consistently by requiring that the 
self-energy of the infinite system and that of the cluster 
be the same. To break the symmetry, frequency indepen- 
dent source fields are allowed on bath sites only. More 
details can be found in . All the calculations are done 
with exact diagonalization at zero temperature (43| . 

For the case of a 2 x 2 plaquette, which we shall con- 
sider throughout this work, the Nambu spinor is defined 
by = (c|p . . . , cj|, ci|, . . . , C4|), and the greek let- 
ters /X, = l,...,Nc label the degrees of freedom within 
the cluster. We compute the cluster propagator Gc by 
solving the cluster impurity Hamiltonian that will be de- 
scribed shorty. Given the Qo on the cluster that results 
from the presence of the bath, we extract the cluster self 
energy from Sc — Qq^ — G~^. Here, 



Gc {t,t') 



Gt(t,t') F{t,t') 
F^{r,T') -Gi(r',r) 



(3) 



is an 8 X 8 matrix, Gf^^^a = -{TCf^a{T)cl^{0)) and Ff^^ = 
— {TCfj_]-{T)c^l (0)) are the imaginary-time ordered normal 
and anomalous Green functions respectively. 

Using the self-consistency condition. 



Go ^(it^n) 



i2n) 



dk G(k, iujn) 



+ ^ciiuJn) (4) 



with 



G(k, zu;„)= iw„ + ^ - £(k) - Ec(ici^„) 



(5) 



we recompute the Weiss field Gq^, obtain the correspond- 
ing bath parameters by minimizing a distance function 
described below, and iterate till convergence. Here i(k) 
is the Fourier transform of the superlattice hopping ma- 
trix with appropriate sign flip between propagators for 
up and down spin and the integral over k is performed 
over the reduced Brillouin zone of the superlattice. 

A 2 X 2 plaquette is embedded in a bath of non- 
interaction electrons. To solve the cluster impurity prob- 
lem, we express it in the form of a Hamiltonian ffimp with 
a discrete number of bath orbitals coupled to the clus- 
ter and use the exact diagonalization technique (Lanczos 



ct ta ct 
f n' fi 



fii^a ni(7 
a 

Here = l,...,Nc label the sites in the cluster and 
Efj,ua represents the hopping and the chemical potential 
within the cluster. The energy levels in the bath are 
grouped into replicas of the cluster {Nc = 4) (two replicas 
in the present case) with the labels m = 1 , ■ • ■ , Nc and 
a = 1,2 such that we have 16 bath energy levels 
coupled to the cluster via the bath-cluster hybridization 
matrix V^^^. Using lattice symmetries we take V^^^ = 
V^Smfi and = e". The quantity A" represents the 
amplitude of superconducting correlations in the bath. 
No static mean-field order parameter acts directly on the 
cluster sites. 

The parameters e", V" and A" are determined by im- 
posing the self-consistency condition in Eq.|4]using a con- 
jugate gradient minimization algorithm with a distance 
function 



(6) 



that emphasizes the lowest frequencies of the Weiss field 
by imposing a sharp cutoff at lUc — 1.5. (Energies are 
given in units of hopping t, and we take H = 1.) The 
distance function in Eq. ^ is computed on the imaginary 
frequency axis (effective inverse temperature, (3 = 50) 
since the Weiss field Go{i(^n) is a smooth function on 
that axis. 

With the bond superconducting order parameter de- 
fined as 



(7) 



we consider d-wave singlet pairing (tp = il)i2 = —'023 — 
■034 — — 04i). The average is taken in the ground state 
of the cluster. 

All quantities depending on wave vector, including self- 
energy, are obtained from the Green function periodiza- 
tion scheme. 

The finite size of the bath in the exact-diagonalization 
technique is an additional approximation to the CDMFT 
scheme. The accuracy of this approximation can be ver- 
ified by comparing the CDMFT solution for the one- 
band Hubbard model with the solution from the Bethe 
ansatz [3]. We have also used this comparison in one 
dimension as a guideline to fix the choice of parameters 
in the distance function {ojc = 1.5 and /3 = 50). These 
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results in one dimension also compare well with those 
obtained using the Hirsch-Fye Quantum Monte Carlo al- 
gorithm as an impurity solver where the bath is not trun- 
cated [45^]. Further, using finite-size scaling for these low 
(but finite) temperature calculations p^, it was shown 
that, at intermediate to strong coupling, a 2 x 2 cluster 
in a bath accounts for more than 95% of the correla- 
tion effect of the infinite size cluster in the single-particle 
spectrum. Because of the finite size of the bath, one also 
needs to use a finite linewidht broadening 77 = 0.125 when 
plotting the figures. 

We can also perform an internal consistency check on 
the effect of the finite bath on the accuracy of the calcu- 
lation. With an infinite bath, convergence insures that 
the density inside the cluster is identical to the density 
computed from the lattice Green function. In practice, 
we find that there can be a difference of ±0.02 between 
the density estimated from the lattice and that estimated 
from the cluster. We display results as a function of clus- 
ter density since benchmarks with the one-dimensional 
Hubbard model show that, with a finite bath and the 
procedure described above, one can reproduce quite ac- 
curately Bethe ansatz results for n (fj.) when the cluster 
density is used. Nevertheless, we should adopt a con- 
servative attitude and keep in mind the error estimate 
mentioned above. 



SUPPLEMENTARY FIGURES AND 
APPENDICES 

Technical comments on Fig. 1 

The off-diagonal self-energy E"„ is extracted from 
the periodized Nambu Green function and plotted 
with 77 = 0.125. The maxima of the plot are marked 
with black dots. The imaginary part of the local spin 
susceptibility x" {^) is calculated on the cluster and 
plotted with the same value of rj. The dots indicate the 
maxima in the limit 77 = 0. The various densities evalu- 
ated on the cluster are plotted for the following values 
of (/i, S) starting from the normal state : (0.25, 0.37) 
(0.375,0.35), (0.5,0.33), (0.625,0.31), (0.75,0.29) 
(0.875,0.26), (1.0,0.24), (1.125,0.22), (1.25,0.20) 
(1.375,0.18), (1.5,0.16), (1.625,0.14), (1.75,0.13) 
(1.875,0.11), (2.0,0.10), (2.125,0.08), (2.25,0.07) 
(2.375,0.05), (2.5,0.04). 



Relation between shift in peak position and 
single-particle gap 

Fig. [3] illustrates how the single-particle gap in the su- 
perconducting state Asc (not necessarily the supercon- 
ducting gap) and the shift d between the position of the 
peaks in S"„ and x" change with doping. 
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FIG. 4: (Color online) As a function of doping, the shift d 
(open circles and blue line) between the first peak in the x" 
and the first peak in E"„ at the antinodal point. Also shown 
as a function of doping, the single particle gap A^c (triangles 
and solid red line) measured from the single-particle density 
of states. On the left panel, U = St, t' /t = 0. On the right 
panel U — 12t and the band parameters are those appropriate 
for YBaaCusOr-i. 



Real part of the anomalous self-energy 

In conventional Migdal-Eliashberg theory, the real part 
of the self energy EJj„ times the quasiparticle renormal- 
ization factor is essentially the gap function. We find 
that this function, illustrated in Fig [5l increases as one 
approaches half-filling, consistent with the increase in the 
single particle gap found earlier [2^ and illustrated in 
Fig. m Sjj„ has weak frequency dependence near zero fre- 
quency only over a range of order J — AP/U for U >8t, 
as can be seen in Fig. [6l If there were a "static" piece to 
the gap, Sj,„ would have a frequency independent compo- 
nent at frequencies larger than J, at least until frequen- 
cies of order U. We find that this is not the case. For the 
t — J model 0, H^l one finds a small instantaneous contri- 
bution to EJin, thus making connection with mean- field 
theories. We show in the main text how mean-field the- 
ories can also be seen as approximations to the present 
approach, even though we do not find an instantaneous 
contribution to pairing. 



The Lehman representation for the Nambu Green func- 
tion allows us to find the following result for the T ~ 
value of Ig (t^) 

Ig H = 1"^) ("^1 Cji |0) e {lo - {E,n - Eo)) 

(8) 

with the Heaviside step function. Excited states |m) 
that contribute have an energy less than oj above the 
ground state |0) . 

For BCS s-wave theory, Ig {oj) can be computed ana- 



7 



o 

'M 



U/t=8 






tvt=o 




- 






— n=0.95 






n=0.89 






n=0.80 






n=0.75 





-2-10 1 2 

(O/t 

FIG. 5: (Color online) Real part of the anomalous self-energy 
"San for U = 8, t' = t" = at the antinodal point. Four dif- 
ferent dopings are presented. Negative contributions appear 
at a frequency of order J nearly independent of doping. 
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FIG. 6: (Color online) Real part of the anomalous self-energy 
Ti'an as a function of frequency u at the antinodal point for 
fixed doping S = 0.11 and different values of i7 = 8t, 12i 
and 16t, t' = t" = represented respectively by solid bleue 
line, short-dashed red line and long-dashed green line. The 
nearly flat part near uj = decreases with J. The range of 
frequencies where is positive also decreases as U increases 
or J decreases. 



lytically. One obtains, using F^^ with i = j, 



sinh"^ (uj/A) - sinh"^ (1) 



sinh ^ (Wc/A) - sinh ^ (1 
X e {lu - A) 9 {uJc - uj) + {uj - LUc) 



(9) 



The results for the d-wave case in the main text were 
obtained by numerical integration and a sharp cutoff. 



FIG. 7: (Color online) Imaginary part of the spin susceptibil- 
ity x" for U = 8t, t' — t" = expressed in cluster momenta 
for underdoping on the top panel {5 = 0.05) and for overdop- 
ing on the bottom panel {S = 0.2). Even in the latter case, 
a sizable (tt, tt) component is left. Recall that the momenta 
refer to averages over a quarter of the Brillouin zone. 



In Ehashbergh theory, the function Ig is 



Ig (lo) = N{0) I Re 
/o 



^{Z (w')w')'-S,„ {cj'f 



duj' 



(10) 

where N (0) is the single-particle density of states at the 
Fermi level, the square root is in the upper half-plane 
and 



Z{lo) = 1 



Sii {uj) -f 1:22 {u}) 
2uj 



ill) 



with Tiii (w) the diagonal components of the self-energy 
in Nambu space, S22 (w) — — Sii (— ^) and San (w) = 
S12 (w). The phase is chosen such that there is no con- 
tribution from the second Pauli matrix in Nambu space 

The anomalous Green function F^j (w) entering the 
calculation of Ig {^) in CDMFT was obtained by Fourier 
transforming the anomalous lattice Green function calcu- 
lated with the band Lanczos approach. We used three 
different values 77 — 0.24,0.18 and 0.12 for the small 
imaginary part that must be added to the real frequency 
to obtain the retarded {uj) . The final result for Ig (w) 
is the extrapolation to 77 ~ 0. This is done to smooth the 
function while preserving as much as possible the asymp- 
totic large frequency value. It differs by only a few per- 
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creases as one approaches hafF-filling. This is an illus- 
tration of the detrimental effect of the pseudogap. The 
self-energy in the dressed Green functions entering the 
calculation leads to a decrease in the number of states 
that can pair near the Fermi level. 



Attractive Hubbard model 



FIG. 8: (Color online) Pairing susceptibility calculated with 
the dressed bubble only, i.e. without vertex corrections. The 
decrease near half-filling illustrates the pair-breaking effect of 
the pseudogap. 



cent from the value of the order parameter calculated on 
the cluster. 



Effect of the self-energy in decreasing the pairing 
tendency as one approaches half-filling 

Fig. [5] shows that in the normal state, the pairing 
susceptibility calculated without vertex corrections de- 



The cutoff frequency enters very clearly in the inte- 
grated off-diagonal spectral weight Eq. Since one 
expects that the attractive (instead of repulsive) Hub- 
bard model should behave more like the BCS model, we 
checked that Ig (lo) calculated with C-DMFT for that 
model docs have the structure of the BCS result for s— 
wave. In other words, it vanishes below the gap, and in- 
creases monotonically until a sharp cutoff frequency that 
depends somewhat on U but is of the order of the band- 
width, as expected from the mean-field solution. There is 
some structure in the frequency dependence that is prob- 
ably caused in part by the finiteness of the bath used in 
the calculation, but does not change the overall trend. 



